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ABSTRACT 



We present a uniform study of the chemical abundances of 12 elements (Si, Ca, Sc, Ti, V, Cr, Mn, Co, Ni, Na, Mg, and Al) 
derived from the spectra of 451 stars observed as part of one of the HARPS GTO planet search programs. Sixty eight of these are 
planet-bearing stars. The main goals of our work are: i) the investigation of possible differences between the abundances of stars 
with and without planets; ii) the study of the possible differences in the abundances of stars in the thin and the thick disc. We 
confirm that there is a systematically higher metallicity in planet host stars, when compared to non planet-hosts, common to all 
studied species. We also found that there is no difference in the galactic chemical evolution trends of the stars with and without 
planets. Stars that harbour planetary companions simply appear to be in the high metallicity tail of the distribution. We also 
confirm that Neptunian and super-Earth class planets may be easier to find at lower metallicities. A statistically significative 
abundance difference between stars of the thin and the thick disc was found for [Fe/H] < 0. However, the populations from the 
thick and the thin disc cannot be clearly separated. 



Key words, stars: abundances - stars: chemically peculiar - stars: fundamental parameters 
disc - galaxy: kinematics and dynamics - galaxy: solar neighbourhood 
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1. Introduction 

The discovery that giant-planet host stars have, on 
average, a high er metallicity than non-planet hosts 
19971) . inspired a series of spectroscopic 



( Gonzalez 



studies for stars studied as part of radial- ve l ocity p lanet' 
hunting pr o gram s (e.g. ISantos et al. 
Laws et al. 2003 ; Fischer fe Valentif 



200C, 120011 L2004; 



2005). Gradually, 



similar abundance studies were made for elements other 
than 



Gonzalez et al. 2001 



Beir ao et al 
20061: 



2005 



20021: Israelian et al. 2003 : Bodaghee et al 



Ecuvillon et al 



Sadakane et al 
20031 



Fischer fe Valentil \200A iGilli et al 
200& iTakedal l2007t ). The same 
abundance difference was observed in these elements, 
suggesting that the probability of finding a giant planet 
is strongly related to the abundance content of its host 
star. 



Interestingly, stars with Neptune and super-Earth- 
class planets might have different trends than the ones 
found for Jovians. There are hints that these planets may 
be easier to find at the typical metallicities of non-planet 
(non-Jovian?) host stars or eyen at lower [Fe/H] values 
(|Udrv et al.ll2006l : ISousa et alJl2008t l. However, there are 
still an insufficient number of planets for compiling a ro- 
bust statistical distribution. 

Some authors also referred to the existence of differ- 
ent [X/Fe] trends in planet-host stars compared to stars 
without planets for the same [Fe/H]. These potential dif- 
ferences were detected fo r some elements (e.g. Si, Ti, Mn 
V, Co, Ni, Na, and Mg -iBodaghee et al.ll2003l; IGilli et al 



* Based on observations collected at the La Silla Parana 
Observatory, ESO (Chile) with the HARPS spectrograph at 
the 3.6m telescope, under the GTO program 072.C-0488. 



20061: iflobinson et aljbood) and are difficult to explain 
with galactic chemical-evolution models. Despite that, the 
results obtained are still contradictory and no definitive 
conclusion has yet been re ached. A review of these trends 
was recently published bv lGonzalez fe Laws! (j2007T ). 

The large amount of spectra used in planet-search sur- 
veys with radial velocity methods may also be employed 
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to study the abundances and kinemat ic properties of the 
stars in the solar neighbourhood (e.g. Santos et al. 2003t 
Fischer fc Valentill2005t iGilli et al.|[2006t lEcuvillon et al 



20071) . This data can be used to study the properties 
of the different populations of the Galaxy and, in par- 
ticular^the differences between the thin and thick disc. 



Fuhrmann et al. 1 19981 ) showed that the thick disc stars 



have a distinctive [Mg/Fe] overabundance compared to 
the stars in the thin disc, for the same [Fe/H], at least 
for metallicities bel ow solar. This fact was confirmed by 
Bensbv et al. (12003) and expan ded for other elements. The 



work of lEcuvillon et alJ (|2007l ). however, did not find any 
clear difference in the element abundance for the different 
populations observed in previous studies. 

In this paper, we derive the abundances of twelve 
species (silica, calcium, titanium, scandium, manganese, 
chromium, vanadium, cobalt, nickel, sodium, magnesium, 
and aluminium), based on spectra of the same HARPS 
guaranteed time observations (GTO) from the planet 
search program. This sample contains 451 stars, of which 
68 are planet hosts and the remain ing 383 stars are d warfs 



with no known orbiting planets (jSousa et alJ 120081 ). We 



also derive the kinematical properties of the sample. The 
main goals of this work are twofold: i) the investigation of 
possible differences between the abundances of stars with 
and without planets; ii) the study of the possible differ- 
ences between the abundances of stars from the thin and 
thick disc of the Galaxy. In Sect. [21 we present the sample 
used in this work, along with the calculation of the galac- 
tic space velocity data and the selection of different groups 
of stars, based on their kinematic properties. The method 
and the selection steps that precede the chemical abun- 
dance determination will be explained in Sect. [3] Section 
@] contains the abundance analysis of the full sample as 
well as the final [X/H] values. It also includes discussion 
of the uncertainties and errors in our methodology as well 
as a comparison of our results with the literature. A dis- 
cussion of our results follows in Sect. where we explore 
the possible differences between stars with and without 
planets in our sample. In Sect. we discuss the differ- 
ences between the abundances of the thin and thick disc 
stars, using the [X/Fe] distributions for [Fe/H] < and 
the [X/Fe] versus [Fe/H] plots. We also use these plots to 
complete a brief study of the trends in the Galactic chem- 
ical evolution. Finally, in Sect. [71 we draw our concluding 
remarks. 

2. The sample 

The spectra for the HARPS "high precision" GTO pro- 
gram were obtained with the HARPS spectrograph at the 
ESO 3.6m telescope (La Silla, Chile). This sub-sample 
consists of 451 stars, selected from the volume-limited 
sample of solar neighbou rhood stars studi ed with the 
CORALIE spectrograph (jUdrv et al.l l2000f) as well as 



a group of planet hosts from the so uthern hemisphere . 
Details about this sample are given in lSousa et al.l (|2008l ). 
We note that we updated the original GTO catalogue 



using data from the extrasolar planets encyclopaedist, 
to take account of th e newly discovered planets orbit- 
ing the stars HD40307 (jMavor et al.ll2008l ). HD47186, and 
HD181433 (Bouchy et al. in prep.). The total number of 
planet-bearing stars in the sample is now 68. Most stars 
in the catalogue are slow rotators, non-evolved, and with 
a low chromospheric activity level. 

The individual spectra were reduced using the HARPS 
pipeline and later combined with IRAF0 after correcting 
for the radial velocities of the stars. The final spectra have 
a resolution of R ~ 110 000 and a signal-to-noise ratio 
ranging from ~ 70 to ~ 2000 depending on the amount 
and quality of the original spectra. Nine tenths of the spec- 
tra have S/N higher than 200 and 50% have a S/N higher 
than 450. 

A solar spectrum was also collected with HARPS using 
solar light reflected by the asteroid Ceres. This spectrum 
has the same spectral resolution as the remaining spectra 
and a S/N ratio of ~ 250. 

The stellar parameters u sed in the present paper were 
determined bv lSousa et alJ (|2008l ). We refer to this paper 
for details. 

2.1. The UVW data and thin/thick disc stars 

The galactic space- velocity components (UVW), anal- 
ysed in Sect. [6l were computed with the trigonometric 
parallax es and pro per motions from the Hipparcos cat- 
alogue (|ESAl Il997h . and by adopting the mean radial 
velocities obtained from the HARPS spectra (courtesy 
of the HARPS GTO team). The Galactic space- velocity 
components were calculated using the procedure from 
Johnson fc Soderbloml (|l987h and corrected for the solar 



motion relative to the local standard of rest (LSR) us- 
ing (Un,Vg,Wn) = (+1 0.00, +5.25, +7.17) km s" 1 from 
Dehnen fe Binnevi (Il998l). T his LSR correction was also 



used bv lBensbv et all (|200l |2005h 

The selection of the thin and thick d isc stars was 
compl e ted us ing the method proposed by iBensbv et al 



(2003 



20051 ). The characteristic velocity dispersions 
(au, &v, &w), the asymmetric drift, V asym , and the ob- 
served fractions of these p opulations in the sol ar neigh- 



bourhood were taken from IBensbv et al. (|2005h . as well 



as the equations to calculate the relative probability of 
a star belonging to a certain population. Each probabil- 
ity (P) is calculated by multiplying its Gaussian proba- 
bility distribution with the observed fraction of the re- 
spective population . Following the recommendations of 
IBensbv et all (|2005l) . if Pthickdisc/Pthindisc > 2, the star 
will be considered a member of the thick disc. Otherwise, 
if Pthickdisc/Pthindisc < 0.6, the star will belong to the 



1 http://exoplanet.eu/ 

2 IRAF is distributed by National Optical Astronomy 
Observatories, operated by the Association of Universities for 
Research in Astronomy, Inc., under contract with the National 
Science Foundation, U.S.A. 
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Table 1. Sample tabic of the parameters used to assign the Galactic population to which each star belongs. Columns 
2, 3, and 4 list the U,V and W velocities relative to the local standard of rest (LSR). Columns 5, 6, and 7 list the 
Gaussian distributions of each population: thin disc, thick disc, and halo, respectively. In Cols. 8 and 9, we provide 
the relative probabilities of the thick disc relative to the thin disc and of the thick disc relative to the halo. The last 
column depicts the probable population where each star belongs. 



Star ID 


Ulsr 


Vlsr 


Wlsr 


fthin 


f thick 


fhalo 


P(thickdisk) 
P(th in disk) 


P(thickdisk) 
P(halo) 


group 






[km s _1 ] 












HD 134060 


26 


-34 


15 


1.75E-06 


5.75E-07 


5.67E-09 


0.04 


6765.43 


thin 


HD134606 


-11 


-16 


1 


5.36E-06 


5.20E-07 


3.80E-09 


0.01 


9134.99 


thin 


HD134664 


27 


2 


-14 


1.92E-06 


2.69E-07 


2.26E-09 


0.02 


7933.64 


thin 


HD134985 


-29 


-65 


-41 


7.14E-09 


2.92E-07 


9.84E-09 


4.55 


1979.36 


thick 


HD134987 


-11 


-35 


28 


7.24E-07 


4.91E-07 


5.62E-09 


0.08 


5824.82 


thin 


HD136352 


-109 


-42 


43 


4.76E-10 


8.89E-08 


4.93E-09 


20.78 


1203.01 


thick 


HD136713 


9 


-13 


-5 


5.15E-06 


4.84E-07 


3.51E-09 


0.01 


9182.53 


thin 


HD136894 


5 


3 





3.76E-06 


3.11E-07 


2.28E-09 


0.01 


9078.65 


thin 
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& 50 
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Fig. 1. Toomre diagram for the full sample (451 stars). 
The stars from the thin disc, transition group, and the 
thick disc are marked by red circles, black triangles, and 
blue squares, respectively. The stars with and without 
planets are shown with full and empty symbols, respec- 
tively. The halo star HD 104006 is represented by a starred 
symbol. 



thin disc. All remaining stars having a probability between 
these values will be included in a transition population. 

The star HD104006 with a Pthickdisc/ Phala = 0.286 
will be considered as belonging to the halo. 

In Fig. [TJ the Toomre Diagram shown represents the 
entire sample. The blue squares, black triangles, and red 
circles correspond to the thick disc, transition and thin 
disc stars, respectively. The stars with and without plan- 



ets are indicated by filled and empty symbols, respectively. 
The halo star HD 104006 is represented by a starred sym- 
bol. A sample of the distributions and probabilities cal- 
culated for each star, as well as the parameters used in 
their calculation, are presented in Table [TJ Among the 
451 stars in our sample, we have 400 stars from the thin 
disc, 29 from the thick disc, and 21 are considered to be 
transition stars that do not belong to any group or have 
a determined origin. According to this adopted criteria, 
only 1 star is from the halo. 



3. Chemical abundances 

The chemical abundances of the studied elements were 
derived using a differential local thermodynamic equilib- 
rium (LTE) analysis. This was c ompleted with the 2002 
version of the MOOGo pr ogram (ISnedenlll973h using the 
'abfind' driver. A grid of iKurucj (|l993l ) ATLAS9 atmo- 
spheres were used as input, along with the equivalent 
widths and the atomic parameters, wavelength (A), exci- 
tation energy of the lower energy level (xz), and oscillator 
strength (loggf) of each line. The atmospheric parame- 
ters, effective temperature (T e ff), surface gravity (log g), 
micr oturbulence ( &), an d metallicity ([Fc/H]) were taken 
from Sousa et al.l ( 2008f ) and derived from the same spec- 
tra used in our study. This provided consistency to our 
work, ensuring that it was more uniform and less affected 
by external errors. The reference abundances used in the 
differ ential analysis are those o f the Sun, and were taken 
from lAnders fc Grevesse (jl989h . 



3.1. The line list 

The initial line list, along with the atomic parameters (in- 
cluding an initial estimate of the oscillator strengths) were 



The source code of MOOG2002 can be downloaded at 



http://verdi.as.utexas.edu/moog.html 
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taken from the VALEn online database (Piskunov et al 



19951 : iKupka et al.l Il999f ). When extracting the lines, we 
used the solar stellar parameters as input: T e ff = 5777 K; 
log g = 4.44 dcx; £ t = 1.0 km/s. We found 2169 lines in 
the spectral region of interest (4500 to 6910 A). After a 
lengthy process of selection and identification of the 'good 
lines', involving the use of t he IRAF "splot" t ool with 
the Kurucz Solar Flux Atlas ( Kurucz et al. 1984f k we ob- 
tained a shorter list of 364 lines. All strongly blended lines, 
as well as lines that were too weak (having Equivalent 
Widths (hereafter EW) below 5mA), or too strong (lines 
with EW > 200 mA) , and those inside the wings of strong- 
lines (e.g. H a (A = 6562.81 A), Hp (A = 4861.34 A), and 
A5172.70, A5183.62 Mg I lines) were excluded in this pro- 
cess. 

The next step involved the automatic measurement 
of the equivalent width of the lines, in a resolution- 
degraded spectrum (R ~110 000) of the Kurucz Solar 
Flux Atlas. This w as completed with the ARES^I program 
(jSousa et al.ll2007l ). Here, we also devised a selection pro- 



Teff S 5300 K 



cedure that involved the comparison of the EWs measured 
in the aforementioned spectrum with those measured in 
the solar reflected light spectrum of the Ceres asteroid, 
obtained with HARPS. Both spectra have a resolution of 
the order of ^110 000. If the relative difference between 
the EWs, for a certain line, measured in the two spec- 
tra, was higher than 10%, the line would be discarded. 
Additionally, we also excluded all lines that ARES could 
not reproduce well due to a poor evaluation of the lines 
in the spectral area of interest, or due to a poor deter- 
mination of the continuum placement. These criteria are 
important because the abundance and especially the os- 
cillator strength determination are very sensitive to varia- 
tions in the EW. The parameters adopted in ARES were: 



smoother 
and rejt = 



4, space = 2, lineresol 



0.996 ffrom ISousa et al 



= 0.07, 
2008). 



minilinc = 5, 
At the end of 



this step, we finished with 284 lines. 

The EW of these lines, along with a model atmosphere 
based on the solar parameters, T e ff = 5777, logg = 4.44, 
logei?e = 7.47, and £ t = 1.0, were used in the calcu- 
lation of the semi-empirical oscillator strengths by us- 
ing an inverse analysis with the MOOG driver "ewfind". 
All input paramete rs for this driver were taken following 



Santos et alj (|2004l ). The new oscillator values are impor- 
tant because the log gf data from VALD might not be 
sufficiently accurate. 

To test the calculated log gf values, we inverted the 
previous analysis using a solar atmosphere model and the 
measured EWs of the solar reflected light spectrum of 
Ceres using the "abfind" driver within MOOG. In general, 
the average abundance values found were within 0.01 dex 
of the expected value, except for Al (A^4 = ±0.04 with 2 
lines), Scl (AA = ±0.03 with 3 lines) and V {AA = ±0.02 
with 14 lines), where A is the abundance of a species. 



Vienna Atomic Line Database 

The Automatic Routine for Line Equivalent Widths is 




S000 5200 5400 5600 5800 6000 6200 6400 6600 6800 
Wavelength (A) 

5300 < Teff < 6000 K 




5000 5200 5400 5600 5800 6000 6200 6400 6600 6S00 
Wavelength (A) 

Teff 1 6000 K 



0.8 
0.6 
0.4 
0.2 


-0.2 
-0.4 



^=|#it= = | ==%f = = =?=! = = =t-t^ 



available at: http://www.astro.up.pt/~sousasag/ares/ 



5000 5200 5400 5600 5800 6000 6200 6400 6600 6800 
Wavelength (A) 



Fig. 2. Process of selection of the stable lines for the neu- 
tral atom of Si, shown as an example. The black crosses on 
the three plots represent the difference between the mean 
abundance of each line (from all the stars) and the mean 
abundance of the element. The error bars represent the 
rms of the mean values. The dashed lines mark the exclu- 
sion value. The open circles correspond to the difference 
values of the cool group (T e ff < 5300 - top plot), the solar 
group (5300 < T eff < 6000 - middle plot), and the hot 
group (T eff > 6000 - bottom plot). 



The final selection was the most important on e. After 
determ ining the abundances for all stars, using the lKurucz 
(|l993l ) model atmospheres calculated with the stellar pa- 
rameters and EWs measured with ARES for each star and 
element, we divided the stars into three distinct groups: a) 
the cool group, containing 152 stars with T e ff < 5300; b) 
the solar group, having 257 stars with 5300 < T e // < 6000; 
and c) the hot group, comprising 42 stars with T e /f > 
6000. This was completed to account for the differences 
in the local conti nuum about each line that are temper- 
ature dependent (jSousa et al.ll2008T ). The ARES parame- 
ters were identical to those of the solar case apart from 
'rejt' which depended on the S/N ratio of each spectra 
according to Sousa et al. ( 20081 ). and 'miniline=2' (ensur- 
ing that all lines measured in the solar spectrum were 
accounted for in the other stars). Then, we calculated the 
difference between the mean abundance of each line (from 
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Table 2. Wavelength (A), excitation potential oscillator strength (loggf) and solar equivalent widths (EW ) 
for the lines selected in the present paper. 



A 


Xi 


logs/ 


EW 


3 


A 




Xl 




logs/ 


EW 


3 


A 


Xl 


logs/ 


EW0 


[A] 


[oV] 




[mA] 


r X i 

[A] 




[cV] 






[mA] 


[A] 


[oV] 




r £ i 

[mA] 


Na I - 2 


lines, 


log e D — 6.33 






5039.96 




0.02 




-1.199 


75 


3 


5737.07 


1.06 


-0.815 


10.4 


5688.22 


2.10 


-0.628 


121 


4 


5064.06 




2.69 




-0.471 


5 


5 


6081.45 


1.05 


-0.692 


14.1 


6154.23 


2.10 


-1.622 


36 


6 


5071.49 




1.46 




-0.797 


28 


8 


6224.51 


0.29 


-1.935 


5.3 


6160.75 


2.10 


-1.363 


54 


3 


5113.44 




1.44 




-0.861 


27 





6251.83 


0.29 


-1.431 


15.0 


Mg 1-3 


lines 


log e = 7.58 






5145.47 




1.46 




-0.622 


37 





6274.66 


0.27 


-1.751 


8.3 


4730.04 


4.35 


-2.234 


69 


6 


5219.70 




0.02 




-2.254 


28 


1 


6285.17 


0.28 


-1.676 


9.5 


5711.09 


4.35 


-1.777 


105 


6 


5490.16 




1.46 




-1.008 


21 


4 


Co I - 9 


lines, lc 


ge„ = 4.92 




6319.24 


5.11 


-2.300 


25 


2 


5503.90 




2.58 




-0.218 


12 


3 


4594.63 


3.63 


-0.279 


12.5 


Al I - 2 


ines. 


logee = 6.47 






5648.57 




2.49 




-0.410 


1(1 


1 


4792.86 


3.25 


-0.080 


32.7 


6696.03 


3.14 


-1.571 


36 


2 


5662.16 




2.32 




-0.123 


23 


5 


4813.48 


3.22 


0.177 


45.9 


6698.67 


3.14 


-1.886 


21 


1 


5739.48 




2.25 




-0.781 


7 


7 


5301.05 


1.71 


-1.950 


19.5 


Si I - 18 


lines 


log e = 7.55 






5766.33 




3.29 




0.326 


9 


6 


5342.71 


4.02 


0.606 


32.3 


5517.54 


5.08 


-2.496 


12 


9 


5965.84 




1.88 




-0.492 


26 


7 


5352.05 


3.58 


0.004 


24.4 


5645.61 


4.93 


-2.068 


35 


8 


5978.55 




1.87 




-0.602 


22 


6 


5359.20 


4.15 


0.040 


9.6 


5684.49 


4.95 


-1.642 


61 


2 


6064.63 




1.05 




-1.941 


8 


3 


5647.24 


2.28 


-1.594 


14.0 


5701.11 


4.93 


-2.034 


37 


7 


6091.18 




2.27 




-0.445 


15 





6814.95 


1.96 


-1.822 


18.8 


5753.64 


5.62 


-1.333 


45 


7 


6126.22 




1.07 




-1.416 


22 


1 


Ni I - 48 


lines, 


og e = 6.25 




5772.15 


5.08 


-1.669 


52 





6258.11 




1.44 




-0.435 


51 


5 


4512.99 


3.71 


-1.467 


19.4 


5797.87 


4.95 


-1.912 


4:! 


9 


6261.10 




1.43 




-0.491 


49 


2 


4811.99 


3.66 


-1.363 


25.6 


5948.54 


5.08 


-1.208 


85 


4 


6599.12 




0.90 




-2.069 


9 


3 


4814.60 


3.60 


-1.670 


16.7 


6125.02 


5.61 


-1.555 


31 


7 


Ti II - 


8 


lines, 


loge 


= 4.99 






4913.98 


3.74 


-0.661 


55.7 


6142.49 


5.62 


-1.520 


33 


1 


4583.41 




1.16 




-2.840 


31 


8 


4946.04 


3.80 


-1.224 


26.2 


6145.02 


5.62 


-1.425 


38 


8 


4636.33 




1.16 




-3.152 


19 


8 


4952.29 


3.61 


-1.261 


32.3 


6195.46 


5.87 


-1.666 


17 


1 


4657.20 




1.24 




-2.379 


49 


3 


4976.33 


1.68 


-3.002 


37.7 


6237.33 


5.61 


-1.116 


61 


1 


4708.67 




1.24 




-2.392 


48 


9 


4995.66 


3.63 


-1.611 


17.9 


6243.82 


5.62 


-1.331 


44 


8 


4911.20 




3.12 




-0.537 


52 


3 


5010.94 


3.63 


-0.901 


48.8 


6244.48 


5.62 


-1.310 


46 


2 


5211.54 




2.59 




-1.490 


32 


8 


5081.11 


3.85 


0.064 


93.5 


6527.21 


5.87 


-1.227 


38 


9 


5381.03 




1.57 




-1.904 


59 


1 


5094.41 


3.83 


-1.108 


30.3 


6721.85 


5.86 


-1.156 


44 





5418.77 




1.58 




-2.104 


49 


5 


5392.33 


4.15 


-1.354 


12.0 


6741.63 


5.98 


-1.625 


15 


5 


Mn I - 


6 


lines. 


log e 


= 5.39 






5435.86 


1.99 


-2.432 


51.7 


Ca I - 12 lines, loge — 6.36 






4502.21 




2.92 




-0.523 


57 





5462.50 


3.85 


-0.880 


41.0 


5261.71 


2.52 


-0.677 


97 


7 


4671.77 




2.89 




-1.567 


14 


8 


5587.87 


1.93 


-2.479 


52.9 


5349.47 


2.71 


-0.581 


94 


6 


4739.11 




2.94 




-0.462 


60 


1 


5589.36 


3.90 


-1.148 


26.7 


5512.98 


2.93 


-0.559 


84 


8 


5377.62 




3.84 




-0.068 


40 


8 


5625.32 


4.09 


-0.731 


37.8 


5867.56 


2.93 


-1.592 


25 


1 


5399.47 




3.85 




-0.104 


38 


5 


5628.35 


4.09 


-1.316 


14.7 


6156.02 


2.52 


-2.497 


9 


6 


5413.67 




3.86 




-0.476 


21 


6 


5638.75 


3.90 


-1.699 


9.8 


6161.29 


2.52 


-1.313 


60 


6 


Cr I - 22 


lines 


log c 


= 5.67 






5641.88 


4.11 


-1.017 


24.1 


6166.44 


2.52 


-1.155 


69 


9 


4575.11 




3.37 




-1.004 


1(1 





5643.08 


4.16 


-1.234 


15.1 


6169.04 


2.52 


-0.800 


92 


2 


4600.75 




1.00 




-1.457 


82 


3 


5694.99 


4.09 


-0.629 


43.1 


6449.82 


2.52 


-0.733 


98 


1 


4626.18 




0.97 




-1.467 


83 


3 


5748.36 


1.68 


-3.279 


28.0 


6455.60 


2.52 


-1.404 


56 


3 


4633.25 




3.13 




-1.215 


10 


4 


5805.22 


4.17 


-0.604 


40.8 


6471.67 


2.53 


-0.825 


91 


2 


4700.61 




2.71 




-1.464 


14 





5847.00 


1.68 


-3.410 


23.0 


6499.65 


2.52 


-0.917 


85 


8 


4708.02 




3.17 




-0.104 


55 





5996.73 


4.24 


-1.010 


20.3 


Sc I - 3 


ines. 


log c = 3.10 






4730.72 




3.08 




-0.345 


46 


3 


6086.29 


4.27 


-0.471 


43.5 


4743.82 


1.45 


0.297 


8 





4767.86 




3.56 




-0.599 


16 


1 


6108.12 


1.68 


-2.512 


65.0 


5520.50 


1.87 


0.562 


6 


6 


4775.14 




3.55 




-1.025 


6 


9 


6111.08 


4.09 


-0.823 


34.2 


5671.82 


1.45 


0.533 


14 


6 


4936.34 




3.11 




-0.343 


45 


6 


6119.76 


4.27 


-1.316 


10.9 


Sc II - 6 


lines 


loge = 3.10 






4964.93 




0.94 




-2.577 


38 


1 


6128.98 


1.68 


-3.368 


25.3 


5526.82 


1.77 


0.140 


70 


3 


5214.14 




3.37 




-0.784 


16 


4 


6130.14 


4.27 


-0.938 


22.1 


5657.88 


1.51 


-0.326 


66 


9 


5238.97 




2.71 




-1.427 


15 


9 


6175.37 


4.09 


-0.534 


49.0 


5667.14 


1.50 


-1.025 


33 


9 


5247.57 




0.96 




-1.618 


82 





6176.82 


4.09 


-0.266 


63.7 


5684.19 


1.51 


-0.946 


37 


2 


5287.18 




3.44 




-0.954 


1(1 


4 


6177.25 


1.83 


-3.538 


14.6 


6245.62 


1.51 


-1.022 


34 


9 


5296.70 




0.98 




-1.373 


93 


1 


6186.72 


4.11 


-0.888 


30.5 


6320.84 


1.50 


-1.863 


8 


1 


5300.75 




0.98 




-2.125 


58 


9 


6204.61 


4.09 


-1.112 


22.0 


Ti I - 30 


lines 


loge = 4.99 






5781.18 




3.32 




-0.886 


15 


5 


6223.99 


4.11 


-0.954 


27.7 


4555.49 


0.85 


-0.575 


63 


7 


5783.07 




3.32 




-0.472 


31 


4 


6230.10 


4.11 


-1.132 


20.6 


4562.63 


0.02 


-2.718 


1(1 


9 


5787.92 




3.32 




-0.183 


46 





6322.17 


4.15 


-1.164 


18.4 


4645.19 


1.73 


-0.666 


22 





6661.08 




4.19 




-0.234 


12 





6327.60 


1.68 


-3.086 


38.6 


4656.47 


0.00 


-1.308 


68 


6 


6882.52 




3.44 




-0.392 


32 


2 


6360.81 


4.17 


-1.145 


18.5 


4675.11 


1.07 


-0.939 


38 


1 


Cr II - 


3 


lines 


log Co 


= 5.67 






6378.26 


4.15 


-0.830 


31.8 


4722.61 


1.05 


-1.433 


18 


7 


4588.20 




4.07 




-0.752 


68 


5 


6598.60 


4.24 


-0.914 


24.9 


4820.41 


1.50 


-0.429 


43 


1 


4592.05 




4.07 




-1.252 


47 


6 


6635.13 


4.42 


-0.779 


23.6 


4913.62 


1.87 


0.068 


511 


2 


4884.61 




3.86 




-2.069 


23 


4 


6767.78 


1.83 


-2.136 


79.2 


4997.10 


0.00 


-2.174 


31 


6 


V I - 7 


lines, log e Q — 


4.00 






6772.32 


3.66 


-0.963 


49.2 


5016.17 


0.85 


-0.657 


63 


2 


5670.85 




1.08 




-0.482 


18 


8 


6842.04 


3.66 


-1.496 


24.2 



all the stars) and the mean abundance calculated using all 
the lines. This was completed for each element, as we can 
see for Si in Fig. [2] The open circles in the top, middle 
and bottom panels, represent data for hot, solar, and cool 
group stars, respectively, corresponding to this difference. 
The mean value of each line is depicted by a black cross 
and the error bars represent its rms. 



We considered all lines with an overall line dispersion 
in any plot of more than a factor of 1.5 higher than the 
rms, to be non-stable and these lines were excluded. The 
dashed lines on the plots of Fig. [5] indicate this exclusion 
value. We also excluded any line with a dispersion in any 
plot clearly beyond the average abundance value. This se- 
lection was done with a careful manual inspection line by 
line and element by element. In a few cases, we chose to 
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relax these rules slightly, especially when there was a line 
close to the exclusion values occurring in only one group. 
Wc also verified each individual line for stray points clearly 
outside of the 2er distribution (due to bad pixels, cosmic 
rays or other unknown effects). These could modify the 
true value of the mean or the dispersion of a certain line 
and lead to a wrong exclusion. After this last selection, 
the final list contained 180 lines, as shown in Table [2 



4. Abundance analysis of the full sample 

As stated before, the chemical abundances of the elements 
were derived by completing a differential LTE analysis rel- 
ative to the Sun. The analysis of the full 451 star sample 
was described in Sect. 13.11 The final abundance for each 
star and element was calculated to be the average value of 
the abundances for all lines detected in a given star and 
element. A set of FORTRAN programs were developed 
to coordinate these calculations and compile an organised 
and easy-to-use database for the output data. A sample of 
our results are presented in Table [3j The complete results 
are available online. 



4.1. Uncertainties and errors 

Uncertainties can introduce different errors in the calcu- 
lated abundances. Random errors can affect the measure- 
ment of EWs for individual lines (e.g. due to noise and 
cosmic-ray hits) and systematic errors can occur, for ex- 
ample, due to blending, NLTE effects, or to a poor choice 
of continuum location. To minimise both types of errors, 
we should use high quality data and as many lines as pos- 
sible for each element. Unfortunately, we were only able 
to select 2 or 3 lines for Na, Mg, Al, Scl and CrII. The 
presence of insufficient data can produce higher internal 
random errors, and might also introduce systematic bias 
that can be extremely difficult to identify in our results. 
Therefore, all conclusions regarding these elements should 
be considered with caution. We further note that we per- 
formed a uniform analysis of the sample, which minimised 
possible errors due to, for example, differences in line lists, 
atmospheric parameters, oscillator strength biases, and 
applied methodologies. 

4.1.1. Random errors 

The typical rms uncertainties for the atmospheric param- 
eters are 24 K for T e ff, 0.04 dex for logg, 0.02 dcx for 
[Fc/H], an d 0.03 km s" 1 for These uncertainties arc 
taken from ISousa et al] ( 20081 ) and characterize the rela- 
tive precision of the derived parameters. 

We tested the sensitivity of the derived abundances to 
the uncertainties in the atmospheric parameters. First, for 
each parameter, we chose three different stars with similar 
atmospheric parameters except for the one being tested. 
Afterwards, we generated new atmospheric models, chang- 
ing the tested parameter by an amount required to assure 




vt= 



[M/H] = 



0.19 
L_ 



-6 -5.5 -5 -4.5 

log (EW/lambda) 

Fig. 3. Example of the calculation of the EP and RW 
slopes for Nil lines for the star HD1461. The squares rep- 
resent the different spectral lines, the long dashed line 
the calculated slope and the short dashed line the average 
value of the abundance (log eps). The values of the slopes 
are -0.010 and 0.020, respectively for the upper and lower 
panel. 

that the sensitivity of the abundances could be clearly ob- 
served (AT eff = ±100 K, A log .g = ±0.3 dex, A& = ±0.5 
km s" 1 , and A[Fe/H] = ±0.3 dcx). The obtained results 
are displayed in Table 2J 

In general, we observe that the neutral elements are 
sensitive to temperature changes, whereas the ions are 
most sensitive to changes in surface gravity. The ions are 
also more sensitive to metallicity changes than the neutral 
elements, although the sensitivity is not as significant as 
the one for either T e ff and log g. 

4.1.2. Testing the stellar parameters 

To verify the validity of the stellar parameters, we tested 
our results in a variety of ways. First, we calculated the 
slopes of the derived abundances of the considered lines as 
a function of both the excitation potential (EP) and the 
logarithm of the reduced equivalent width (RW) of the 
Nil lines. In this way, we verified whether the excitation 
equilibrium enforced on the Fel lines of every star and 



used to obtain the stellar parameters ( Sousa et al. 2008h 



was applicable to other species. In other words, we expect 
the slopes to be as small as possible. 

We can see an example of this (star HD1461) in Fig. [31 
Wc have chosen nickel because its lines cover a wide range 
of both EP and RW, thus allowing a reliable determination 
of the slopes. 

In Fig. IU we plot the slopes of EP, obtained for every 
star, as a function of the stellar parameters (T e ff, logg, 
£t, and [Fc/H]). From the analysis of the plots, we observe 
no discernible trends except in the case of the T e ff plot, 
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Table 3. Sample table of the derived abundances of the elements, rms and number of measured lines (n) for each star. 



Star ID 


[Si/H] 


rms 


n 


[Ca/H] 


rms 


n 


[Scl/H] 


rms 


n 


[ScII/H] 


rms 


n 


HD12617 


0.15 


0.10 


17 


0.05 


0.14 


10 


0.48 


0.31 


3 


0.14 


0.11 


6 


HD 13060 


0.06 


0.07 


16 


0.03 


0.06 


12 


0.12 


0.11 


3 


-0.01 


0.05 


6 


HD 13724 


0.24 


0.02 


17 


0.22 


0.02 


11 


0.30 


0.02 


3 


0.31 


0.07 


6 


HD 13789 


-0.04 


0.08 


16 


0.01 


0.12 


11 


0.31 


0.25 


3 


-0.13 


0.13 


6 


HD 13808 


-0.14 


0.07 


16 


-0.20 


0.11 


11 


-0.05 


0.11 


3 


-0.22 


0.08 


6 


HD 14374 


-0.02 


0.02 


16 


-0.03 


0.04 


11 


0.00 


0.05 


3 


-0.07 


0.07 


6 


HD 14635 


0.02 


0.10 


18 


0.05 


0.14 


11 


0.49 


0.30 


3 


-0.04 


0.11 


6 


HD 14680 


-0.16 


0.05 


17 


-0.12 


0.11 


12 


-0.01 


0.12 


3 


-0.21 


0.07 


6 


HD 14744 


-0.11 


0.04 


18 


-0.13 


0.06 


11 


0.01 


0.17 


3 


-0.15 


0.08 


6 


HD 14747 


-0.24 


0.01 


17 


-0.23 


0.01 


11 


-0.22 


0.03 


3 


-0.22 


0.04 


6 



Table 4. Abundance sensitivities of the studied elements to changes in the stellar parameters (temperature with 
AT e ff = ±100 K, surface gravity with A logy = ±0.3 dex, metallicity with A[Fc/H] = ±0.3 dex, and microturbulence 
with A£ t = ±0.5 km s" 1 ). 



Star ID 


Teff [K] 


log g 


[Fe/H] 


ft [km s" 1 ] Na Mg Al 


Si 


Ca 


Scl 


Sell 












Temperature with AT eJf = ±100 K 












LLU 1UOU / 


L ± 1 4:0 


4 49 


-0.02 


0.90 ±0.10 ±0.02 ±0.07 


T0.06 


-i-n 1 9 


-i-n 1 4 


zcfl 09 




Hr>9flfil fi 
nijzuuiu 


OOliJ 


4 4 1 } 


0.01 


0.94 ±0.06 ±0.06 ±0.05 


to.oi 


3ZU.UO 


3ZU. 1U 






HD20945 


6118 


4.50 


0.03 


1.21 ±0.06 ±0.06 ±0.05 


±0.03 


±0.07 


±0.09 


±0.01 










Surface gravity with Alogp = ±0.3 dex 












HD19556 


5676 


4.03 


0.06 


1.11 T0.05 tO.05 tO.02 


±0.01 


T0.05 


T0.01 


±0.13 




HD11720 


5667 


4.32 


0.22 


1.01 T0.07 tO.06 tO.02 


±0.01 


T0.07 


to.oi 


±0.12 




HD20260 


5658 


4.49 


0.18 


1.02 =p0.07 tO.06 tO.02 


±0.01 


T0.08 


to.oi 


±0.12 










Microturbulence with A£ t = ±0.5 km s 1 












HD9828 


5381 


4.42 


-0.26 


0.64 =p0.02 tO.03 t0.02 


to.oi 


=F0.07 


T0.02 


=F0.06 




HD18956 


5726 


4.41 


-0.24 


0.95 T0.02 =p0.04 =p0.01 


to.oi 


T0.08 


to.oi 


T0.07 




HD4444 


5999 


4.37 


-0.22 


1.26 T0.02 tO.03 tO.OI 
Metallicity with A [Fe/H] = ±0.3 dex 


T0.02 


T0.08 


to.oi 


T0.09 




HD7874 


5778 


4.46 


-0.67 


1.03 T0.03 ±0.00 T0.01 


±0.01 


±0.00 


to.oi 


±0.06 




HD9670 


5845 


4.39 


-0.18 


1.04 ±0.00 ±0.01 =F0.01 


±0.02 


±0.00 


to.oi 


±0.08 




HD20220 


5757 


4.47 


0.29 


1.01 ±0.04 ±0.03 ±0.01 


±0.05 


±0.03 


±0.01 


±0.11 




Star ID 


Teff [K] 


log g 


[Fe/H] 


ft [km s^ 1 ] Til Till V 


CrI 


CrII 


Mn 


Co 


Ni 










Temperature with AT eff = ±100 K 












HD10567 


4748 


4.42 


-0.02 


0.90 ±0.14 =p0.02 ±0.15 


±0.09 


T0.08 


±0.06 


±0.02 


T0.02 


HD20616 


5519 


4.43 


0.01 


0.94 ±0.12 =p0.01 ±0.12 


±0.08 


T0.04 


±0.07 


±0.06 


±0.05 


HD20945 


6118 


4.50 


0.03 


1.21 ±0.09 ±0.01 ±0.10 


±0.08 


T0.02 


±0.07 


±0.07 


±0.06 










Surface Gravity with Alogp = ±0.3 dex 










HD19556 


5676 


4.03 


0.06 


1.11 =p0.01 ±0.12 =p0.01 


T0.02 


±0.11 


T0.02 


±0.01 


±0.01 


HD11720 


5667 


4.32 


0.22 


1.01 T0.01 ±0.12 =p0.01 


T0.03 


±0.11 


T0.05 


±0.03 


±0.01 


HD20260 


5658 


4.49 


0.18 


1.02 =p0.02 ±0.12 =F0.01 
Microturbulence with Aft = ±0.5 km s" 


T0.04 
-l 


±0.11 


T0.04 


±0.02 


±0.02 


HD9828 


5381 


4.42 


-0.26 


0.64 =p0.08 tO.06 t0.03 


T0.06 


T0.07 


T0.07 


T0.05 


T0.05 


HD18956 


5726 


4.41 


-0.24 


0.95 tO.03 tO.08 tO.01 


T0.06 


T0.09 


T0.04 


T0.03 


T0.03 


HD4444 


5999 


4.37 


-0.22 


1.26 tO.03 tO.08 tO.01 
Metallicity with A [Fe/H] = ±0.3 dex 


T0.04 


T0.12 


T0.04 


T0.03 


T0.03 


HD7874 


5778 


4.46 


-0.67 


1.03 T0.01 ±0.06 T0.01 


to.oi 


±0.04 


to.oi 


±0.00 


±0.00 


HD9670 


5845 


4.39 


-0.18 


1.04 =p0.01 ±0.07 =F0.01 


to.oi 


±0.06 


±0.00 


±0.01 


±0.02 


HD20220 


5757 


4.47 


0.29 


1.01 ±0.00 ±0.11 ±0.01 


±0.02 


±0.08 


±0.03 


±0.03 


±0.05 



where we observe that cooler stars (T e ff < 5000 K) have 
a systematic bias away from the expected values. A higher 
dispersion is also seen, on average, for cooler stars. This 
might be due to the stronger line blending in cooler stars 



as well as to the fact that the computed logy/ values 
(derived from the solar spectrum) may be inadequate for 
these stars. The same effect was noticed in the RW plots. 
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Fig. 4. EP slopes as a function of T e //, logp, and [Fe/H] 
for Ni. 
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Fig. 5. Plots of [CrI/Crll] as a function of T e ff, logg, £ t 
and [Fe/H]. 



We also plotted the [CrI/Crll] value as a function of 
the stellar parameters to ensure t hat the ionisation equi- 
librium enforced on the Fell lines (jSousa et al. 2008 ) was 
acceptable to other elements. This is depicted in Fig. [5] 
We expect that [CrI/Crll] should be independent of these 
parameters. 

As seen in Fig. [5j there is a divergence in the expected 
ratio for stars with T e ff < 4900 K. We are unable to 
discern any significant trend in the plots of the other pa- 
rameters. This result (and that found for Ni) raises some 
doubt about the abundance values derived for stars with 
lower temperatures (T e ff < 4900 K). 

Finally, we plotted the abundance ratios [X/Fe] as a 
function of the stellar parameters. We did not detect any 
significant trends except for the T e ff plot, which is pre- 
sented in Fig.H The slopes of [X/Fe] with T eff per 1000 
K arc listed in Table El 

Wc observe significant trends with T e f f for Til, Scl, V, 
Co, and Al. It is clearly evident that the trends drift away 
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Fig. 6. [X/Fe] vs. T e ff plots. The dots represent the stars 
of the sample. The solid lines depict the linear fits of the 
data. The vertical purple dotted lines indicate the cutoff 
temperature when appropriate. 



from the expected values, except in the cases of Co and 
Al, are heavily influenced by the cooler stars for which the 
abundance might have been overestimated due to blend- 
ing effects; it may also be sensitive to either deviations 
from excitation or ionisation equilibrium, or to problems 
associated with the differential analysis. We must indeed 
emphasise that the oscillator strengths were calculated for 
the Sun. To account for these effects, we decided to es- 
tablish a cutoff temperature: T cuto ff = 5000 K for Til, 
T cutoff = 5100 K for Scl, and T cutoff = 5300 K for V. 
Only stars of temperature above these values are consid- 



Table 5. Slopes of [X/Fe] ratio as a function of T e // per 
1000K 



Species 


Slope(T)±rms 


Species 


Slope(T)± rms 


Si 


-0.004±0.051 


Crl 


-0.053±0.026 


Til 


-0.192±0.081 


CrII 


0.039±0.048 


Till 


0.062±0.076 


Co 


-0.128±0.067 


Scl 


-0.269±0.128 


Ni 


-0.012±0.034 


Sell 


0.050±0.065 


Na 


-0.017±0.066 


Ca 


-0.045±0.069 


Mg 


-0.009±0.089 


Mn 


-0.034±0.093 


Al 


-0.161±0.080 


V 


-0.452±0.086 
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ered in the remainder of this paper, for the aforementioned 
elements. We note this when we analyse our results. 

4.2. Comparison of the abundances with the literature 

To test the reliability of our results, we co mpare the de 



rived a bundances with those ob t ained by Bcns by et al 



(l20f)5h IValenti fc Fischer! (I2005T ). lOilli et all (j2006h . and 
Takedal (l2007ri 7 for stars in common, as seen in Fig.[7J In 



this way, it was possible to verify qualitatively any sys- 
tematic errors in our abundance determinations. 





Til ^y^ 

y /// ^^ 


Till y . 


/ Tefl>5100 


< ' ' 1 ' ' ' ' 1 ' ' ' ' 1 ' ' ' y 


y ' y 

c.a y 
' ' 1 ' ' ' ' 1 ' ' ' ' 1 ' ' ' y 


1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 V 

Mn S 
Inn, 


CrI y 


on y 


/ Teff>5300 


^mImmImmIimi,, 




"i ' ' ' 1 ' ' ' ' 1 ' ' ' ' 1 ' ' \ y 


Mr y 





Fig. 7. Comparison of t he derived abu n dance s with those 
listed by other authors : Bensbv et al.l ( 2005! ) (red trian- 
gles'). I Valenti fc Fischer (120051) ( black crosses), iGilli et al.l 
jgOod) (blue dots) and [Takedal (|2007h (green stars). The 



elements are identified in the top left corner of each plot. 



It is evident that, in general, our results agree with the 
[X/H] obtained by other authors. However, w e observe a 



syste matic und erabundance trend in C a, CrI (jGilli et al 
|2006[ ). and Ti (|Valenti fc Fischer! 120051) and a systematic 
overabundance trends in Mn dCilli et al.l 12006! ) . Mg, and 
Al (jGilli et all [20061 ; iBensbv et all 12005! ) . We emphasise 
the existence of a pop ulation of stars with a systematic 
underabundance in Si (j Valenti fc Fischer! l2005h and veri- 
fied that this is not a planet-host population. We do not 
know the origin of these differences. However, the differ- 
ences are, in general, small and systematic, implying that, 
from a relative point of view, all studies agree well. 



5. Abundances in planet-hosts 

Although we study a relatively small number of planet- 
host stars (68), this number is sufficient to observe 
whether there are any discernible differences in the abun- 
dances of stars with and without planets. 

We note that five planet-host stars (HD4308, 
HD27894, HD111232, HD114729, and HD330075) are con- 
sidered to belong to the thick disc. This means that ap- 
proximately 7% of the planet-host stars in our sample be- 
long to the thick disc and that approximately 17% of the 
thick disc stars arc planet hosts. These values must be re- 
garded as lower limits, since (some?, most?) planets may 
still be awaiting confirmation in the entire sample. It is 
important to note that there was no bias in the choice of 
our sample towards or against thick-disc planet-host stars. 

5.1. The [X/H] Distributions 

The [X/H] distributions of planet and non-planet hosts 
are depicted in Fig. [5] The stars with and without planets 
as well as the stars only with Neptunes and super- Earth^ 
are represented by a dashed blue, solid red, and dotted 
black line, respectively. The vertical lines on the top of 
each histogram represent their average value. 

We observe a clear metallicity excess for giant planet- 
host stars in all elements we study. These results are 
in good agreement with previous sim i lar studies (e.g . 
Bodaghee etaTll2003l iBeirao et alj |2005l IGilli et al]|2006l: 



Bond et al.M2006t iTakedaf 20071 ). Our results are also con 
sistent with pr evious studies of iron abundances (e.g 



Gonzaled Il998 



Santos et al. 2001, 2003 



Gonzalez et al 
2004 



2001 



2005 



Laws et all 12003 



Fischer fc Valenti 



20051 ). as expected. Interestingly, we also find that among 
thick disc members, stars with planets are more metal- 
rich (<[Fc/H]> = -0.16) than "single" stars (<[Fc/H]> 
= -0.33). 

Table [6] lists the average values of [X/H] along with 
their rms, the number of stars used in their determina- 
tion, and the difference of averages between stars with 
and without planets. These differences range from 0.17 for 
Al to 0.27 for Mn. The last four columns list the average 
[X/H] ratio of Neptunian hosts and hosts with Neptunes 
and super-Earths only, along with the number of stars 
used in their determination. 

The difference in the average abundances of stars with 
and without giant plane ts is, on aver age, slightly above 
the values measured by ( Takeda 20071) by 0.03 dex. T his 
difference is also above the results of Bond et al. I (|2006l) for 
al l elements. Our results are, on average, simila r to those 
of iBodaghee et al.l (|2003l ) and lGilli et all (|2006! ). The dis- 
persion in this differences for all elements is the smallest 
of all authors, which probably attests to the quality of our 
results. 



6 HD4308, HD40307, HD69830, HD47186, HD160691 and 
HD181433. The first three stars host neptunians and super- 



earth type planets only. 
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[X/H] [X/H] [X/H] 



Fig. 8. [X/H] distribution of the different elements. The stars with and without planets are represented by a dashed 
blue, and solid red line respectively. The stars that exclusively host neptunians and super-earth planets are represented 
by a black dotted line. The mean value is pictured by a dashed (with planets), solid (without planets), and dotted 
(neptunian and super-earth hosts only) vertical lines on the top of each distribution. The neptunian and super-earth 
distribution was set smaller for clarity. 

It is interesting to see that in most histograms (except host stars are asymmetrical: there is an increase in [X/H] 
for Mg) the distributions of the abundances in planet- to a maximum value and shortly afterwards there is a 
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Table 6. Average abundance values for stars with and without planets, along with their rms, the number of stars used 
in their determination, and the difference between the two groups. The last four columns list the average abundance 
values of neptunian hosts and hosts with neptunes and super-earths only, along with the number of stars used in their 
determination. 



Species 


Planet hosts 




Non-planet hosts 


Difference 


Neptunian 


hosts 


Neptunian 


hosts only 


(X) 


([X/H]) 


rms 


N 


([X/H]) 


rms 


N 


of Averages 


([X/H]) 


N 


([X/H]) 


N 


Si 


0.13 


0.17 


68 


-0.08 


0.21 


383 


0.21 


0.08 


6 


-0.16 


3 


Til 


0.16 


0.14 


62 


-0.01 


0.18 


321 


0.17 


0.12 


4 


-0.05 


2 


Till 


0.11 


0.15 


68 


-0.09 


0.20 


383 


0.20 


0.04 


6 


-0.18 


3 


Scl 


0.20 


0.19 


61 


-0.01 


0.22 


288 


0.21 


0.14 


4 


-0.07 


2 


Sell 


0.14 


0.19 


68 


-0.09 


0.23 


383 


0.23 


0.12 


6 


-0.17 


3 


Ca 


0.10 


0.14 


68 


-0.08 


0.18 


383 


0.18 


0.05 


6 


-0.14 


3 


Mn 


0.11 


0.24 


68 


-0.15 


0.31 


383 


0.27 


0.05 


6 


-0.28 


3 


CrI 


0.11 


0.17 


68 


-0.10 


0.22 


383 


0.21 


0.06 


6 


-0.18 


3 


Crll 


0.07 


0.17 


68 


-0.14 


0.22 


383 


0.21 


-0.01 


6 


-0.26 


3 


V 


0.18 


0.20 


56 


-0.02 


0.25 


243 


0.20 


0.15 


4 


-0.05 


2 


Co 


0.17 


0.21 


68 


-0.05 


0.25 


383 


0.22 


0.19 


6 


-0.12 


3 


Ni 


0.12 


0.20 


68 


-0.12 


0.24 


383 


0.23 


0.07 


6 


-0.22 


3 


Na 


0.11 


0.20 


68 


-0.11 


0.23 


383 


0.22 


0.10 


6 


-0.20 


3 


Mg 


0.13 


0.16 


68 


-0.07 


0.20 


383 


0.19 


0.12 


6 


-0.08 


3 


Al 


0.15 


0.17 


68 


-0.03 


0.20 


383 


0.17 


0.17 


6 


-0.04 


3 


Fe* 


0.10 


0.17 


68 


-0.12 


0.23 


383 


0.22 


0.03 


6 


-0.24 


3 



* [Fe/H] values from Sousa et al. (2008). 



significant cutoff in distribution. This cutoff is located at 
[X/H] ~ 0.3 for Ca, and Crll, [X/H] - 0.4 for Si, Til, 
Till, Scl, Sell, CrI, Ni, Na, Mg, and Al and [X/H] - 0.5 
for Mn, V, and Co. The cutoff might suggest that this is 
the metallicity lim it of the solar neighbourhood stars (e.g. 



Santos et alj 120031 ). since most of the planet hosts are at 



the high metallicity end of the sample. 

Some planet-host distributions appear to be slightly 
bimodal (Mn, V, and Na). We cannot explain this effect 
which might be due only to small number statistics. 

Regarding stars with only Neptunian planets (N=3), 
we observe that they have the lowest [X/H] mean value 
for each element. We note that there are only 3 stars 
with Neptunes, and 3 stars that are Neptunian and Jovian 
hosts. Therefore, any results regarding the Neptunian and 
super-Earth planet hosts must be seen as preliminary. 
Despite that, the values for the Neptunian and super- 
Earth-only hosts confirm that the probability of find- 
ing one of these planets may be even higher in metal- 
poor stars, and they appear to exhibit a different trend 
from the Jovian-host stars. Our results for s tars hosting 
only neptun ian-li ke planets agree wit h those o f lUdrv et al 
( 20061 120071) and lSousa et~aTT (|2008l) for [Fe/H]. 



Figure [9] illustrates the [X/H] distribution in a differ- 
ent perspective: the histograms of the number of stars with 
planets relative to the total number of stars of each bin 
(0.1 dex). We observe that there is a general increase in the 
percentage of stars with planets for all elements, with in- 
creasing [X/H], after a low metallicity plateau. This obser- 
vation ag rees with the findings of o ther authors such as, for 
examp le. lSantos et all (|200lLl2Q04h . and lFischer fe Valent " 
( 20051) for [Fe/H]. 




Fig. 9. Percentage of stars with planets as a function of 
[Fe/H]. The element label and the used cutoff temperature 
(when required) are located at the upper left corner of each 
plot. 
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Fig. 10. [X/Fc] vs. [Fc/H] plots for the elements in study. The blue triangles and the red circles represent the stars 
with and without planets, respectively. The intersection of the dashed lines indicate the solar value. Each element is 
identified in the upper right corner of the respective plot and the cutoff temperature is indicated in the lower right 
corner when appropriate. 



5.2. The [X/Fe] versus the [Fe/H] plots: stars with 
and without planets 

The [X/Fe] versus [Fe/H] plots are traditionally used to 



identify its different stellar populations (thin disc, thick 
disc , and halo), which have distinct abun dance ratios (see 
e.g. iBensbv et al. I l2003t iFuhrmaiml l2004t ) . In this section, 
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Table 7. Average abundance values for the stars from the thin disc, transition and thick disc populations with [Fc/H] 
< 0, along with their rms. The last two columns depict the difference of the average abundance and the KS probabilities 
between the thin and the thick disc populations. 



Species 


Thin Disc 


Intermediate 


Thick Disc 


Difference 


KS 


(X) 


([X/Fe]) 


rms 


([X/Fe]} 


rms 


<[X/Fe]> 


rms 


of Averages 


probability 


Ol 


0.04 


0.04 


0.10 


0.07 


0.14 


0.08 


n in 
u. ru 


9 Qo no. 


Til 


0.10 


0.07 


0.19 


0.10 


0.22 


0.11 


0.12 


6.7e-06 


Till 


0.03 


0.07 


0.13 


0.10 


0.19 


0.09 


0.16 


2.7e-ll 


Scl 


0.07 


0.06 


0.14 


0.11 


0.15 


0.07 


0.08 


1.5e-06 


Sell 


0.02 


0.06 


0.11 


0.07 


0.14 


0.07 


0.12 


2.1e-ll 


Ca 


0.06 


0.06 


0.12 


0.12 


0.14 


0.07 


0.08 


8.9e-07 


Mn 


-0.05 


0.07 


-0.12 


0.12 


-0.17 


0.08 


0.12 


5.6e-10 


CrI 


0.03 


0.04 


0.04 


0.06 


0.03 


0.03 


0.00 


6.5e-01 


CrII 


-0.01 


0.05 


0.01 


0.06 


0.04 


0.10 


0.05 


2.7e-05 


V 


0.05 


0.06 


0.10 


0.07 


0.11 


0.08 


0.06 


3.1e-03 


Co 


0.05 


0.08 


0.11 


0.09 


0.10 


0.06 


0.05 


8.3e-05 


Ni 


-0.01 


0.03 


0.02 


0.02 


0.02 


0.03 


0.03 


2.6e-04 


Na 


-0.00 


0.06 


0.06 


0.09 


0.06 


0.06 


0.06 


6.1e-07 


Mg 


0.05 


0.08 


0.17 


0.10 


0.21 


0.11 


0.16 


2.1e-09 


Al 


0.10 


0.10 


0.21 


0.14 


0.22 


0.09 


0.12 


8.8e-06 



we analyse them only to determine whether there are any 
differences in the abundances of stars with and without 
planets, for the same value of [Fc/H]. 

Figure [TOl shows the [X/Fe] values for the two samples 
as a function of [Fc/H]. The blue triangles and the red 
circles represent stars with and without planets, respec- 
tively. The red dashed lines target the solar value. The 
cutoff temperature adopted is indicated where appropri- 
ate (see Sect. l4.1.2|) . 

The plots do not show any clear differences between 
the stars with and without planets for a given [Fe/H] 
value. We also completed a Kolmogorov-Smirnov (KV) 
test of the [X/Fc] distributions for [Fe/H] > without 
identifying any significant differences between the planet- 
host and the non-planct-host populations. The results 
are in general agre ement with previous si milar studies 



(2001) 



such a s, for instance. I Gonzalez etal. I (l200lf).lTakeda et al 



Beirao et al 



Sadakane et al.l 
(l2005h. 



20021). 



Bodaghee et al 
20051) ~ 



(l2003h 



To avoid systematic errors, we also completed the anal- 
ysis only for stars with T e ff = T Q ±300 K. These plots 
provided us with a more accurate picture of the trends 
and differences among the two groups of stars: in a differ- 
ential analysis such as this, the closer the star temperature 
is to the solar temperature, the more accurate, on average, 
the derived abundance is. A KS test was also completed 
for the [X/Fe] distribution of this temperature-restricted 
population. We did not find any new trend. 

Some outliers may however exist in our sample. In all 
cases, the abnormal abundances may be explained by their 
high dispersion in individual line abundances, by a high 
value of EP and/or RW slopes (see Sect. 14. lH?)) . or as being 
due to non-LTE systematic effects for the highest temper- 
ature stars. 

We did not identify any abundance anomalies in the 
stars previously described in the literat ure: the under- 
abund ance of Ca in HD209100 derived bv lBodaghee et al" 



.Fischer fc Valentil „ 
(120061) . and lTakedal (|2007f) 7 which encourages confidence 
in our results. 

The reported cases of potenti al trends described b 



Gilli et all <|2003h and the overabundance of Mg and Al in th e 



host star HDT687 46 de scribed by 



Laws et al.l (|2003l ). and I Gilli et al 



Sadakane et~al1 (|2002h . 
( 20061 ) was not found. 



other authors such as, for i nstance. iSadakane et al.l (|200 



for vanadium and c obalt, Bodaghee et al. ( 2003f ) 



, for Ti, 

Mn, V, and Co, an dlGilli et all (120061) for V, Co, Mg, and 



Al, were not found. LRobinson et alJ ( 20061 ) claimed to have 
found cl ear and significant overab undances of [Si/Fe] and 



[Ni/Fe]. iGonzalez fc Laws! (|2007t ) also found that there 



were strong indications that Al and Si abundances are sys- 
tematically smaller for the planet-host stars in the higher 
metallicity region, and that the Ti abundance exhibits the 
opposite trend, hinting that the abundances of Na, Mg, Sc 
and Ni might have some differences between planet-hosts 
and non-planet-hosts. We could not identify any of these 
trends clearly in our data. 



6. Thin and thick disc 

In this section we perform a detailed analysis of the [X/Fe] 
distributions of three different populations: thin disc, thick 
disc, and the intermediate population, for [Fc/H] below 
the solar value. Above [Fc/H] > 0, we were unable to 
identify any clear difference between thin and thick disc 
stars. The chosen criteria for selecting stars in each popu- 
lation and the necessary data for determining membership 
were described in Sect. 12.11 We must note that we only 
have 6 thick disc stars with [Fe/H] above the solar value. 
HD27894 is the thick disc member with the highest metal 
content of approximately 0.2 dex. 
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6.1. The [X/Fe] distributions 

The [X/Fe] distributions for [Fe/H] < are depicted in 
Fig. [TlJ Each element and ion is identified at the top left 
corner of the plots, along with the adopted cutoff tem- 
perature when relevant. The solid red, dotted black, and 
dashed blue lines represent the thin, intermediate, and 
thick disc populations, respectively. The mean value of 
each population is depicted by the vertical lines above the 
distributions. 

Table [7] lists the average values of [X/Fe] for the three 
populations and their rms values. The last two columns 
depict the difference of the average abundance and the KS 
probabilities between the thin and thick disc populations. 

From the analysis of the histograms and the associ- 
ated table, we can observe that in the metallicity domain 
[Fc/H]<0, thick disc stars exhibit higher average [X/Fc] 
values for most of the species. However, for CrI, there is no 
difference between the two populations and for Co and Ni, 
this difference is small (0.03 and 0.04 dex respectively). We 
can also see that for all elements there is at least a partial 
superposition between the thin and the thick disc popula- 
tions. The manganese histograms also display a different 
trend from the remaining elements: the values of thin and 
thick disc populations invert, and, this time, the thick disc 
has an underabundance relative to the thin disc. The dif- 
ference between both populations in this element is 0.12 
dex. 

The K-S probabilities imply that all the stars from the 
thick and the thin disc are indeed from separate groups, 
with the exception of CrI, which has a 65% probability of 
belonging to the same population. We note that the re- 
sults based on Crll data differ from those of CrI. However, 
CrII only exhibits two or three lines and, therefore, its re- 
sults might be affected by more significant bias. The mean 
values of the transition population are situated between 
those of the other populations, but its average oscillates 
somehow for some elements. 

It is difficult to conclude anything about the origin of 
this intermediate population. This is also valid for any 
statistical feature of the transition or the thick disc pop- 
ulations that might appear in the histograms: the data 
set is too small in size to enable reliable statistical mea- 
surements to be made. Nevertheless, we can observe that 
most thick disc histograms appear to be bimodal (Si, Til, 
Till, Ca, Mn, V, Co, Na, Mg, and Al): this may indicate 
that there are two populations in the thick disc with dis- 
tinct [X/Fe] values. Alternatively, the observed behaviour 
could be due to the fact that some stars from the thin 
disc may be misclassificd as thick disc stars. The stars be- 
longing to the lower [X/Fe] branch of the thick disc pop- 
ulation have typically lower values of Pthickdisc / Pt hindisc 
when compared to the upper-branch thick-disc stars (for 
Mn this tendency is opposite, as expected). 

With these results we confirm that the thin and thick 
disc populations are chemically and statistically distinct 
(at least in the region below solar metallicity). Besides the 
differences already discovered by iFuhrmann et aL (1998) 



for Mg and bv lBensbv etaD (|2003l ) for Si, Ca, Ti and Al, 
we also found that Sc, Mn, Co and Na display the same 
separatio n between the aforem entioned populations. The 
results of Bensbv et al. ( 2005f) for Na differ from ours: for 
this element, they found that the thin-disc stars might 
be more abundant (have higher [X/Fe]) than thick-disc 
stars. We found the opposite trend. Regardi ng Ni, we can 
only s ay that there is a hint, also found bv iBensbv et al 



(]2003l ). that the thick-disc stars might be more abundant 
than thin-disc stars at [Fe/H] < 0. 



6.2. The [X/Fe] versus [Fe/H] plots: thin and thick 
disc populations 

The [X/Fc] versus [Fe/H] plots are depicted in Figs. [12] 
and [131 The blue, black crossed, and red empty circles 
represent the thick-disc, transition, and thin-disc stars, 
respectively. The dashed lines indicate the solar value. In 
Fig- [ISl only stars with T e ff±300K are shown: this allows 
us to obtain a clearer picture of the trends in the two 
groups. 

Following the results of the [X/Fe] distributions, in Fig. 
[TT1 we verify that, in the [X/Fe] versus [Fe/H] plots, the 
abundance difference between the thin and thick disc pop- 
ulations is not clearly defined. We can observe that, for 
[Fc/H] < 0, the thin and the thick disc populations are 
mixed, and for a fixed value of [Fc/H], wc cannot con- 
clude that the populations are fully separated. Instead, 
for the a elements (Si, Ca, Ti, and Sc) and magnesium, 
the plots reveal a bifurcation into two different groups of 
stars, containing members from the three considered pop- 
ulations, that converges at [Fe/H] between -0.2 dex and 
the solar value. There are some hints that this might also 
be the case for V, Co, Na, and Al. These trends are clearer 
in Fig. [13] The slope in the branch with the higher [X/H] 
is always steeper than the slope of the lower branch. This 
means that the change in the [X/H] ratio of the higher 
branch is faster for [Fc/H] below the solar value. 

We did not identify any differences between the planet 
and non planet-hosts among the thick disc members for 
a given [Fe/H] value. However, we found some indication 
that metal-poor ([Fe/H] ^ —0.2) planct-host stars origi- 
nate prefe r ential ly in the thick disc, as also suggested by 
Havwoodl ( 20081 ). Three of the five planet hosts in this 
metallicity domain are classified as thick disc stars, one as 
a transition object and the remainder as member of the 
thin disc. However, we note that we have only five stars 
in these conditions, and therefore, are unable to confirm 
reliably the validity of this result. 

IBensbv et all (|2003l) and IBensbv et all (|2005f ) con- 
cluded that there is a clear separation between stars from 
the thick and thin disc populations. No reference to the 
existence of the bifurcation discussed here was given by 
these authors. We are still unable to ascertain its origin, 
although as discussed in the previous section, it may indi- 
cate that some thin-disc members have been misclassificd. 
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6.3. The [X/Fe] versus [Fe/H] plots: galactic chemical 
evolution 
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thin disc stars, respectively. The dashed lines mark the solar value. Each element is identified in the upper right corner 
of the respective plot and the cutoff temperature is indicated in the lower left corner when appropriate. 



generati ons of stars with the exception of the lightest el- 
ements (|McWilliamlll997l ). Here, we use the [X/Fc] versus 
[Fc/H] relations, in Figs. [T2l and [T3l to study the chemical 
evolution of the galaxy, assuming that [Fe/H] can be used 



as a time variable. It is difficult to say whether the heavy 
scatter found in some plots are cosmic (astrophysical) or 
due to errors. 
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In these figures, we consider the stars with and 
without planets as one group, but we identify stars in 
the thick and thin-disc, as well as in an intermediate 
group (see Sect. 12. ip . The analysis is completed ele- 
ment by element and compared with other studies of 
chemical abundances between stars with and without 



lanet s, namely, Sadakane et al. I2OO2I). Bodaghee et al 



.20031). iBeirao et alj (I2005D. iFischer fc Valentil (120051). 
Gilli et all (|2006l) . iGonzalez fc Laws! (|2007t) . and iTakeda 
d2007|) (hereafter SAD, BOD, BEI, FV05, GIL, GZ07, 
and TAK, respec t ively) , as well as other reference stud 



ies bv lChen et all (|200of ). lAllende Prieto et"al] (|2004 ) , and 
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Bensbv et all (|2005h (hereafter CHOO, AP04, and BEN, 
respectively) . 

6.4. The alpha elements 

It is believed that the alpha elements (Si, Ca, Ti, and Sc 
in this study) are produced mostly in massive stars that 
explode as Type II supernova (SNe II), although Type I 
supe rnova (SNe la) might als o yield some of these elements 



(e.g. iThielemann et al.ll2002fl 



We can observe that for all elements in this group, 
there is a bifurcation in the value of [Fc/H] below solar. 
The upper branch descends from overabundances of ~ 0.3 
and ~ 0.4 dex toward solar abundances at [Fe/H] around 
the solar value, while the lower branch has a much slower 
descent. The abundance then remains at solar value and 
at some point it ascends, descends, or remains constant 
for [Fc/H] above the solar value, depending on the ele- 
ment: this might be observational evidence that the a el- 
ements are produc ed in different amounts by different SN 
(jMcWilliaml Il997l ). We should note, however, that tita- 
nium is overabundant relative to the other a elements at 
[Fe/H] below the solar value, and that its slope (in this 
metallicity region) is steeper than for calcium and sili- 
con. The Sc abundances, for [Fe/H] < 0, seem to slowly 
decrease or remain constant until [Fe/H] reaches ~ —0.3 
dex, where it has a downturn toward solar values. 

6.4.1. Silicon 

We can see that there exists a clear bifurcation in a region 
of [Fe/H] ranging from -0.85 to 0.0 dex. This can be seen 
more clearly in Fig. [T3] The [Si/Fe] value of the upper 
branch decreases until [Fe/H] ~ —0.2 dex, before joining 
with the lower branch, that has a much gentler descent. 
Then, the [Si/Fc] ratio remains constant at solar value, 
with a slight upwardly trend starting at [Fc/H] ~ +0.2. 

The results of other authors are, in general, similar to 
our own. However, the plateau described also in CHOO ap- 
pears to start much earlier than in our case, and in BEN 
there seems to be no plateau. BOD, FV05, and GZ07 
obtained slightly lower abundances (~ 0.05 dex) in the 
plateau region. The upward trend at [Fc/H] > is not 
evident in SAD, FV05, GIL, or TAK. The AP04 measure- 
ments are systematically higher than ours ( by more than 
~ +0.1 dex) and the upward slope is far more significant 
than in our results. No author refers to the bifurcation 
observed here. 



6.4.2. Calcium 

Calcium shows an uniform downward trend for the entire 
metallicity range. It begins from abundance values of ~ 
+0.3 at [Fc/H] ~ -0.85 dex, and decreases to [Ca/Fe] 
values of ~ —0.1 at [Fe/H] ~ +0.4, equalling the solar 
value at solar [Fe/H]. The branches of the bifurcation are 
more difficult to observe and they seem to join at [Fe/H] ~ 



—0.3 dex. A slight plateau might exist for solar abundance 
values in the metallicity region between —0.2 and the solar 
value. This is more evident in Fig. 1131 

The trends obtained for calcium are similar to most 
studies found in the literature, such as CHOO, BOD, GIL, 
TAK and GZ07. However, the abundances obtained by 
SAD, AL4, and BEN do not continue to descend for 
[Fe/H]> 0. The results of GIL and GZ07 present the same 
trend observed here, but the [Ca/Fe] values are systemat- 
ically lower than ours. No author refers to the presence of 
a bifurcation. 

6.4.3. Titanium 

The upper branch of both Til and Till plots decrease 
until [Fc/H] reaches ~ 0.0-0.2 dex, where it joins the 
lower branch. It then settles there for the remainder of 
the [Fe/H] range, at solar [Ti/Fe] value. Since the plots of 
Ti exhibit considerable scatter, it is easier to observe the 
trends in Fig. [13] We can see that the higher branch pop- 
ulation has a [Ti/Fe] value between ~ +0.35 and ~ 0.0 
dex, in a region of [Fe/H] ranging from —0.85 to ~ 0.0 
dex. 

A similar behaviour in the titanium trends is found 
in the studies of CHOO, SAD, BOD, GIL, BEN, TAK, 
and GZ07. The plateau region in BOD, however, appears 
to be slightly below the solar value. The observed trend 
in FV05 has a continuously decreasing trend throughout 
the [Fe/H] range. The data from AP04 exhibit an upturn 
in the [Ti/Fe] values above solar metallicity. We did not 
detect the weak upturn in Till reported by TAK. Again, 
no author refer to the observed bifurcation. 



6.4.4. Scandium 

There is considerable scatter in the points in Fig. [12] for 
the Sc plot. The trend is easier to observe in Fig. [13] In the 
Scl panel, we observe a slow decrease in the abundance of 
both branches until they join at [Fe/H] ~ —0.1. It remains 
constant at the solar value in the region of —0.1 < [Fc/H] 
< +0.1, and there is then an increase in the abundance 
for [Fc/H] > +0.1 dex. The Sell plot is similar to that 
of Scl, but there appears to be no plateau: the [ScII/Fc] 
value decreases until the [Fe/H] reaches the solar value 
and then begins to increase toward the limit of the [Fe/H] 
range. 

The results of BOD, GIL, and GZ07 differ slightly from 
ours. The results are similar for [Fe/H] < 0, but then we 
observe an increase in the abundance with [Fc/H], which 
they do not observe. However, BOD suggests that there 
might be an increase for iron-rich stars, but this is incon- 
clusive due to scattering. This increase was also noticed by 
GIL and GZ07, although it is elusive. Our figures suggest 
a much stronger trend. The results of SAD and TAK are 
similar to ours, although the Scl and Sell results of the 
latter author do not show the upturn or show a very weak 
one for metallicity values above the solar value. AP04 also 
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reports the same trends, but their results have a system- which some authors propose is caused by NLTE effects 
atic overabundance when compared to ours. No branch is (e.g. BOD, GIL), 
found in the referred works. 



6.5. The iron peak elements 



6.5.3. Cobalt 



It is believed that the iro n peak elements are mos tly cre- 
ated in SNe la explosions ( Thielemann et aT1l2002l ). In our 
plots, we can observe that CrI and Ni (for [Fe/H] < 0) fol- 
low iron in lockstep, which implies that they should have 
the same origin. On the other hand, vanadium and cobalt 
have trends more similar to the a elements. Manganese 
seems to mirror the behaviour of the alpha elements. 



6.5.1. Chromium 

In the CrI plot, the chromium abundances are constant 
around the solar value, as expected by the theory. The 
scatter is small when compared to the other elements. 
We note that the CrII plot has a weak, continuous down- 
ward trend but we only detect 2 or 3 lines per star, which 
means that the abundance determination is more affected 
by blending effects or to some other unknown systematic 
error. 

The CrI results of CH00, SAD, BOD, BEN, and TAK 
are similar to ours. However, we observe a slight system- 
atic underabundance in the results of BOD. The results 
of GIL for CrI suggest a downward trend similar to that 
in our CrII plot and their values have a systematic un- 
derabundance when compared with ours. Regarding our 
CrII plot, we find that it is different from the CrI plots of 
every author. The only work, besides ours, with a CrII de- 
termination is TAK: it shows a systematic upward trend 
throughout the metallicity range, in disagreement with 
our results. 



6.5.2. Vanadium 

Vanadium has a behaviour similar to the alpha elements, 
which might indicate a common origin. However, we do 
not see a clear bifurcation, although there are some hints 
that it might exist (see Fig. [T3|) . The [V/Fe] values slowly 
decline toward solar metallicity up to [Fe/H] ~ —0.1. The 
trend then becomes inverted and we can observe an in- 
crease toward the upper end of the [Fe/H] region. We can 
easily verify that there is a lot of scatter. We should note 
that the vanadium abundances could not have been deter- 
mined if we had not removed the stars with T e tf > 5300 
K (sec Sect. l4.1?2|) . Therefore, it is impossible to reach any 
important conclusion about the true trend or the origin of 
this particular clement. 

Our results for vanadium differ from those of CH00 
and SAD, who reported that the vanadium abundances 
reflect that of Cr and Ni. However, the results of BOD, 
GIL and, TAK are similar to ours. Nevertheless, we note 
that almost every author reports high levels of scattering, 



Cobalt has a large scatter throughout the entire metal- 
licity range. The abundance value decreases slowly until 
it reaches [Fe/H] ~ —0.1 dex. The [Co/Fe] ratio then re- 
mains constant around the solar value until [Fc/H] reaches 
the solar metallicity. Toward higher metallicities, there is 
a slow but clear increasing trend. In Fig [Ljl the abun- 
dance seems to be constant all the way to the solar value 
and a bifurcation similar to that found for the a elements 
might also be present. Cobalt has trends similar to both 
vanadium and the alpha elements. This clement might be 
subjected to NLTE effects (see Sect. 14.1 .2"]) . Cobalt, as well 
as vanadium, are odd-Z nuclei but, despite that, they have 
trends similar to those of the a elements. 

The general trends observed by SAD, BOD, GIL and 
TAK are present in our plots. AP04 also exhibit a similar 
result. However, the upturn observed by AP04 for high 
[Fe/H] is steeper than those observed here. All abundance 
determinations, including ours, have considerable scatter. 
It is thus difficult to draw any strong conclusions from this 
analysis. 



6.5.4. Nickel 

The trends seen in Figs. [T^] and Q]|] are all similar for 
this element. The [Ni/Fe] value remains roughly at so- 
lar value until [Fe/H] ~ —0.3 dex. It then has a small 
decrease, almost step-like, and remains in a plateau un- 
til [Fe/H] reaches the solar value. From this point on, it 
has two strong upturns at solar [Fe/H] and at ~ +0.15 
dex. Between these values and from ~ +0.15 to the edge 
of the metallicity range, we can observe the existence of 
two small plateau with solar [Ni/Fe] ratio and [Ni/Fe] 
~ +0.05, respectively. Nickel abundances have a low dis- 
persion. We speculate that each one of these latter 'steps' 
might correspond to events that initiated an increase in 
the production of Galactic Ni. 

Our results are similar to those of BEN, although a lit- 
tle different in their fine details: they only report one up- 
turn at solar metallicity We note an approximate agree- 
ment with the results of CH00, SAD, FV05, and TAK, 
since they all present evidence of the same general trends 
(slow decrease up to the solar metallicity, followed by a 
weak upturn). Although BOD, GIL, and GZ07 report sim- 
ilar trends, their results have an abundance slightly below 
ours in the entire [Fc/H] range (~ —0.05 dex). The same 
trends are also observed by AP04 but their abundance 
values are higher than most results and their slopes are 
also steeper. 
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6.5.5. Manganese 

In contrast to all other elements, we observe a uniformly 
increasing trend for [Mn/Fe] values from ~ —0.4 to 0.1 dex 
throughout the entire [Fe/H] range. In Figs. fT2l and IT3| we 
have the impression that two different branches exist for 
[Fe/H] values below ~ 0.1 dex. The manganese undcr- 
abundance in this region seems to mirror the a-elemcnt 
overabundance . 

There is good agreement between our results and those 
of SAD and BOD. However, they do not find any hint of 
the existence of two separate branches. The GIL and TAK 
results differ from ours: they show a plateau for [Fe/H] 
below the solar value. 



~ —0.1 dex. The abundance has an increasing trend for 
[Fc/H] > +0.2. There is considerable scatter in Fig. [T3J 
It is easier to follow the trends in Fig. [TSJ There is also 
a slight hint in Fig. Q1J] that there might be a bifurcation 
below solar [Fe/H]. The trend of Na is similar to the alpha- 
element trends, but the overabundance values at negative 
[Fe/H] are not so pronounced. 

A similar behaviour was found by SAD as well as by 
BEI, FV05, GIL, BEN, TAK, and GZ07. CH00 reported 
a constant abundance value throughout the entire [Fe/H] 
range which disagrees with our results. No author men- 
tions the existence of a plateau in the intermediate [Fe/H] 
region. No bifurcation is reported in any paper. 



6.6. The light metals: Mg, Na and Al 

Sodium and aluminium arc understood to be 



mostly in the cores of massive stars (jChen et al 



c reate d 
2000h . 



that later become SNe II. We can observe that the trends 
of Na and Al are similar to those of the alpha elements: 
we might even consider them to be mild alpha elements, 
even if we take into account the fact t hat they have an odd 
number of protons ( McWilliamlll997j ). It is still difficult to 
differentiate the sources of Na and Al. Magnesium is pro- 
posed to be created only in SNe II, and can be considered 
a member of the alpha elements. 

6.6.1. Magnesium 

The [Mg/Fe] branches exhibit a decreasing trend until 
[Fc/H] ~ —0.1, although the upper branch has a much 
steeper descent. The branches then join and the abun- 
dance remains constant at solar value. We observe con- 
siderable scatter in the Mg plots, due to the fact that we 
used only 2 or 3 lines per star in the determination of 
its abundance. Magnesium is an interesting element: it is 
predicted that it c an only be created in SNe II explosions 
( Chen et al. I l2000h . but we can see that there is no con- 
tinuo us decline in magnesium as we would expect in this 
case (jChen et al.l 2000: B cnsbv et al.l 120031) . Therefore, a 



contribution of SNe la or other unknown factors may be 
causing the enrichment of Mg in the solar neighbourhood. 

The obtained trends arc similar to those measured by 
SAD, AP04, BEI, BEN, GIL, and GZ07, but none of 
these authors refers to a bifurcation below solar [Fc/H] 
value. CH00 detect a similar trend but their plateau starts 
much earlier, at —0.3 dex, and the abundance value of this 
plateau for the entire range of [Fc/H] is higher than ours 
by 0.1 dex. TAK, on the other hand, have a systematic un- 
derabundance in the same plateau above solar metallicity 
when compared with our results. 

6.6.2. Sodium 

Sodium abundance values slowly decrease until they reach 
the solar value at [Fe/H] ~ -0.2 dex. They then remain 
constant for a metallicity range from [Fc/H] ~ —0.2 to 



6.6.3. Aluminium 

The trends for aluminium are very similar to those of the 
alpha elements. One might say, at least from an observa- 
tional point of view, that aluminium is an alpha clement. 
Its abundance decreases until [Fe/H] ~ -0.2 dex. It then 
remains at constant abundance (^+0.05 dex) until [Fc/H] 
reaches +0.1 dex. Afterwards, there is an upturn followed 
by a plateau with [Al/Fe] ~ +0.05 in the [Fe/H] region 
between 0.2 and 0.4 dex. The trends are easier to see in 
Fig. [131 where we can also see the presence of a possible 
bifurcation similar to the one in the Mg plot. Again, we 
observe more significant scatter than for sodium. 

Our results are similar to those obtained by BEI, GIL, 
and BEN. CH00 and SAD obtained different results from 
ours: the abundance remains constant at solar value and 
for CH00, it undergoes an increase toward higher metal- 
licities at which it reaches ^—0.3 dex. 



7. Concluding remarks 

We derived the abundances of 12 elements (Si, Ca, Sc, 
Ti, V, Cr, Mn, Co, Ni, Na, Mg, and Al) in a sample of 
451 field dwarfs from the HARPS "high precision" GTO 
planet search program, of which 68 planets are known 
to harbour planetary companions. The spect roscopic pa- 



ramet ers used in this study were derived by ISousa et al 

(|20oah . 



We compared our results with those presented in other 
works to ensure consistency and reliability in our analy- 
sis. The [X/H] distributions as well as the plots of the 
abundance ratios [X/Fe] versus [Fe/H] were analysed in 
comparing the trends of stars with and without planets. 

We derived the galactic space- velocity components and 
used them to identify kincmatically the origin of the stars 
of our sampl e , acc ordin g to the method de veloped by 



Bensbv et ail (|2003l ) and iBensby et all (|2005h . We found 



that the sample consists of 400 stars from the thin disc, 
29 stars from the thick disc, and 21 from a transition pop- 
ulation that exhibits characteristics in common with both 
populations. One star is from the halo. 

We confirmed that an overabundance in giant-planet 
host-stars is clear for all the studied elements. This may 
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agree with previous suggestions that the efficiency of plan- 
etary formation correlates strongly with the metal con- 
tent of the host star. In other words, the probability of 
finding a star with giant planets is higher in stars with 
high [X/H] ratios. The plot of the percentage of stars with 
planets as a function of the [X/H] ratio also provides us 
with this information. This dependence is also predicted 
by studies of core-accr etion modelling (see e.g. Ilda fc Lin 
2004l:lBenz et al.ll2006l) . We also confirmed that stars host- 
ing only neptunian-like planets may be easier to detect 
around stars with similar or even lower metallicity than 
non-planet hosts. This suggests that lower mass planets 
may have a different kind of formation than Jovian plan- 
ets, which is also predicted by the aforementioned core- 
accretion studies of planetary formation. 

We do not find any evidence of differences between 
stars with and without planets for the same [Fc/H] in the 
plots of [X/Fc] versus [Fe/H]. This is in general agreement 
with previous studies. The stars that harbour planetary 
companions simply appear to be in the high metallicity 
tail of the distribution, following the same trends as the 
stars without planets. This result favours a primordial ori- 
gin as the preferred cause of the observed high mctallicities 
in planet-host stars. 

We found a clear distinction between thick disc and 
thin disc populations, for [Fe/H] < 0, when we anal- 
ysed the [X/Fe] distributions. This extends to all ele- 
ments except chromium and nick el, and confirms and ex - 
tends the general res ults f ound bylFuhrmann et alJ (|1998T ) , 
Bensbv et all (|2003h . and lBensbv et all (|2005l )~ 

Interestingly, however, we observed that in the plots of 
[X/Fc] versus [Fe/H] the populations of the thin and thick 
disc arc mixed for the same [Fc/H] throughout the entire 
metallicity range. We were unable t o identify any sepa ra- 
tion between the po pulations as in iFuhrmarml ( 2004 ) or 
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Instead, the plots of the alpha el- 
ements and Mg reveal, for [Fe/H] < 0, a bifurcation of 
two different groups of stars, both with stars from the 
three populations, that converge at [Fe/H] close to the so- 
lar value. This bifurcation might also exist for V, Co, Na, 
and Al, but it is unclear due to scatter. In some plots (Si, 
Ca, Sc, Ti, Mg and Al), the two populations of this bif- 
furcation seem to have different evolutionary dynamics: 
the one with lower abundances (dominated by thin-disc 
stars) has a shallow slope (if any), while the higher abun- 
dance population (dominated by thick-disc objects) has a 
far steeper slope, implying a more dynamic evolution. 

In general, the trends that we found agree with the 
literature references, with the exception of the aforemen- 
tioned bifurcation. Nevertheless, we note that the tradi- 
tional separations between alpha elements, iron peak ele- 
ments, and others might need to be reformulated as the el- 
ements tend to have different trends, at least in the metal- 
licity range of our study. Our results suggest that the ori- 
gin of these elements is more complex than predicted by 
theory. Some elements exhibit a considerable scatter (V, 
Co, Na, Al), although we were unable to ascertain whether 



this scatter is cosmic in origin or due to errors in the anal- 
ysis. 
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